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We study the dynamics of phase transitions out of equihbrium in weakly coupled 
Qh! scalar field theories. We consider the case in which there is a rapid supercooling 
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from an initial symmetric phase in thermal equilibrium at temperature Tj > Tc to 
a final state at low temperature Tj ~ 0. In particular we study the formation and 
growth of correlated domains out of equilibrium. It is shown that the dynamics of the 
process of domain formation and growth (spinodal decomposition) cannot be studied 
in perturbation theory, and a non-perturbative self-consistent Hartree approximation 
is used to study the long time evolution. We find in weakly coupled theories that the 



size of domains grow at long times as £,D{i) ~ ^t^(O). The size of the domains and 
the amplitude of the fluctuations grow up to a maximum time ts which in weakly 
coupled theories is estimated to be 
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with ^(0) the zero temperature correlation length. For very weakly coupled the- 
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ories, their final size is several times the zero temperature correlation length. For 
strongly coupled theories the final size of the domains is comparable to the zero 
temperature correlation length and the transition proceeds faster. 



11.10.-z;11.9.+t;64.90.+b 



Typeset Using REVTEX 



2 



I. INTRODUCTION AND MOTIVATION: 



Phase transitions play a fundamental role in our understanding of the interplay between 
cosmology and particle physics in extreme environments. It is widely accepted that many 
different phase transitions took place in the early universe at different energy (temperature) 
scales and with remarkable consequences at low temperatures and energies. In particular 
broken symmetries, and possibly the observed baryon asymmetry in the universe. 

Phase transitions are an essential ingredient in inflationary models of the early universe 
The importance of the description of phase transitions in extreme environments was 
recognized long time ago and efforts were devoted to their description in relativistic quantum 
field theory at finite temperature . For a very thorough account of phase transitions in 
the early universe see the reviews by Brandenberger Kolb and Turner |]1T| and Linde 

The methods used to study the equilibrium properties of phase transitions are by now 
well understood and widely used, in particular field theory at finite temperature and effective 
potentials ||14 |. 

These methods, however, are restricted to a static description of the consequences of 
the phase transition, but can hardly be used to understand the dynamics of the processes 
involved during the phase transition. In particular, for example, the effective potential 
that is widely used to determine the nature of a phase transition and static quantities like 
critical temperatures etc, is irrelevant for the description of the dynamics. The effective 
potential corresponds to (minus) the equilibrium free energy density as a function of the 
order parameter. This is a static quantity, calculated in equilibrium, and in particular to 
one loop order it is complex within the region of homogeneous field configurations in which 
V"{(j)) < 0, where V{(j)) is the classical potential. This was already recognized in the early 
treatments by Dolan and Jackiw p]. 

In statistical mechanics, this region is referred to as the "spinodal" and corresponds to 
a sequence of states which are thermodynamically unstable. 

At zero temperature, the imaginary part of the effective potential has been identified 



with the decay rate of this particular unstable state ||26|| . 

The use of the static effective potential to describe the dynamics of phase transitions 
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has been criticized by many authors, among them Mazenko, Unruh and Wald |T^. These 
authors argued that phase transitions in typical theories will occur via the formation and 
growth of correlated domains inside which the field will relax to the value of the minimum 
of the equilibrium free energy fairly quickly. It is now believed that this picture may be 
correct for strongly coupled theories but is not accurate for weak couplings. 

The mechanism that is responsible for a typical second order phase transition from a 
initially symmetric high temperature state is fairly well known. When the temperature 
becomes lower than the critical temperature, long-wavelength fluctuations become unstable 
and begin to grow and the field becomes correlated inside "domains" , the order parameter 
(the expectation value of the volume average of the field), remaining zero all throughout the 
transition. 

Recently, de Vega and one of the authors, have studied the influence of these instabilities 



in the evolution of the order parameter out of equilibrium [16] in the case when a non-zero 
(but small) initial value of the order parameter was assumed. 

In this work we continue the study of the non- equilibrium aspects of second order phase 
transitions in typical scalar field theories, with a view towards a deeper understanding of the 
dynamics of phase transitions in inflationary scenarios of the early universe. In particular 
trying to describe the process of domain formation and growth in the case in which the initial 
state is symmetric and in equilibrium at a temperature higher than the critical temperature 
and cooled down below the critical temperature. 

Although there have been several attempts to study the time evolution of the scalar field 
either in flat spacetime or de Sitter space [p!7| -|21||, to our knowledge there has not as yet 
been a consistent treatment of the dynamics of domain formation and growth. 



Recently, Kolb and Wang ||22| have reported on an equilibrium study of the static prop- 
erties of domains produced in late-time phase transitions, but it becomes a pressing issue to 
understand the process of domain formation and growth, especially the time scales involved 
and the size of the correlations. 

Although our initial motivation, and ultimate goal, is to study the dynamics of the phase 
transition in an expanding universe, in this article we report our studies on the dynamics 
of the phase transitions in Minkowski space. We do not attempt in this article to study 
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the non-equilibrium properties in the necessarily more complicated setting of inflationary 
cosmologies, and restrict ourselves to introducing the methods of non-equilibrium quantum 
statistical mechanics and apply them to the study of formation and growth of domains in 
flat spacetime. 

We would like to point out at this stage, that the situation under consideration is very 
different from the classical description of the process of spinodal decomposition in statistical 
mechanics. The classical approach to spinodal decomposition is based on a "coarse grained" 
time dependent Landau- Ginzburg equation, which is flrst order in time and purely dissipa- 
tive. Thermal fluctuations are usually introduced as a Langevin noise, typically uncorrelated 
23| , |2^ , that obeys the fluctuation dissipation relation. 

In our case we are studying a quantum field theory, the Heisenberg fleld equations are 
second order, non-dissipative (in Minkowski space), and both quantum and thermal fluctu- 
ations are present in the initial state (density matrix), furthermore, as is typical in these 
scalar theories, the order parameter is not conserved. 

This article is organized as follows: in section 2 we present our arguments that suggest 
that phase transitions in expanding cosmologies must be studied away from equilibrium for 
weakly coupled theories. We emphasize that the important long-wavelength modes that 
become unstable below the critical temperature and whose dynamics is relevant for the 
process of phase separation and domain growth will easily be out of equilibrium during the 
transition for weakly coupled theories. 

In section 3, we introduce our model and the methods of non-equilibrium statistical 
mechanics as applied to the description of the dynamics of the phase transition. 

In section 4, we analyze the real time correlation functions in zeroth order and obtain 
the flrst quantitative expressions, a scaling law for the size of the domains and the growth of 
the amplitude of the fluctuations. In section 5, we carry out a perturbative calculation and 
show quantitatively that the dynamics cannot be studied within a perturbative framework. 

In section 6 we introduce a non-perturbative self consistent Hartree approximation to 
study the evolution of correlations and growth of domains. We provide an analytic and 
numerical analysis of the process of domain growth and establish a scaling law for the size 
of the domains at long times, and an estimate for the maximum size of the domains for very 
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weakly coupled theories. 

We summarize our findings and pose further questions in section 7. 

II. THE CASE FOR A NON-EQUILIBRIUM DESCRIPTION: 

Before entering into the technical details, let us sketch the arguments that suggest that 
when the temperature is very near the critical temperature, the relevant dynamics must be 
studied out of equilibrium. 

As in any dynamical process, in order to try to describe the time evolution of the system, 
one must first try to determine the typical time scales involved in the different dynamical 
processes. This understanding becomes more pressing when one tries to understand the 
dynamics of phase transitions. 

In particular, is it possible to describe the phase transition in an environment in which the 
temperature is changing at some particular rate, in local thermodynamic equilibrium?. To 
address this issue one must compare the typical coUisional relaxational rates of the particles, 
to the rate of change of the temperature. 

The typical coUisional relaxation rate for a process at energy E in the heat bath is 
given by r{E) — r~'^{E) n{E)a{E)v{E) with n{E) the number density of particles with 
this energy E, a{E) the scattering cross section at this energy and v{E) the velocity of 
the incident beam of particles. The lowest order (Born approximation) scattering cross 
section in a typical scalar theory with quartic interaction is (t{E) ^ X^/E^. At very high 
temperatures, T » m$, with m$ the mass of the field, n{T) fa T^, the internal energy 
density is U/V ^ T^, and the average energy per particle is {E) f» T, and v{T) fs 1. 
Thus the typical coUisional relaxation rate is r(T) ^ A^T. In an expanding universe, the 
conditions for local equilibrium will prevail provided that r(T) 3> {a{t)/a{t)) = H{t), with 
a{t) the FRW scale factor. If this is the case, the collisions occur very quickly compared to 
the expansion rate, and particles will equilibrate. This argument applies to the coUisional 
relaxation of high frequency (short- wavelength) modes for which k ^ m$. This may be 
understood as follows. Each "external leg" in the scattering process considered, carries 
typical momentum and energy k,E T > T^. But in these typical theories ~ m$/\/A- 
Thus for weakly coupled theories Tc ^ m$. 
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To obtain an order of magnitude estimate, we concentrate near the phase transition at 
T ^ T^^ lO^^Gev, H ^ T^/Mpi ^ lO^^T, this imphes that for A > 0.01 the conditions 
for local equilibrium may prevail. However, in weakly coupled inflaton models of inflation, 
phenomenologically the coupling is bound by the spectrum of density fluctuations to be 
A ~ 10"^^ — 10"^"^ []T^,|TB[. Thus, in these weakly coupled theories the conditions for local 
equilibrium of high energy modes may not be achieved. One may, however, assume that 
although the scalar field is weakly self-coupled, it has strong coupling to the heat bath 
(presumably other fields in the theory) and thus remains in local thermodynamic equilibrium. 

This argument, however, applies to the collisional relaxation of short wavelength modes. 
We observe, however, that these type of arguments are not valid for the dynamics of the 
long-wavelength modes at temperatures below Tc, for the following reason. 

At very high temperatures, and in local equilibrium, the system is in the disordered 
phase with ($) = and short ranged correlations, as measured by the equal time correlation 
function (properly subtracted) 



As the temperature drops near the critical temperature, and below, the phase tran- 
sition occurs. The onset of the phase transition is characterized by the instabilities of 
long-wavelength fluctuations, and the ensuing growth of correlations. The field begins to 
correlate over larger distances, and correlated domains will form and grow. If the initial 
value of the order parameter is zero, it will remain zero throughout the transition. This 
is the process of spinodal decomposition, or phase separation. This growth of correlations, 
cannot be described as a process in local thermodynamic equilibrium. 

These instabilities are manifest in the equilibrium free energy in the form of imaginary 
parts, and the equilibrium free energy is not a relevant quantity to study the dynamics. 

These long wavelength modes whose instabilities trigger the phase transitions have very 
slow dynamics. This is the phenomenon of critical slowing down that is observed exper- 
imentally in binary mixtures and numerically in typical simulations of phase transitions. 
The long wavelength fluctuations correspond to coherent collective behavior in which de- 
grees of freedom become correlated over large distances. These collective long-wavelength 



($(f,t)$(0,t)) ^T2exp[-|f|/e(T)] 





(2.2) 
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modes have extremely slow relaxation near the phase transition, and they do not have many 
available low energy decay channels. Certainly through the phase transition, high frequency, 
short wavelength modes may still remain in local equilibrium by the arguments presented 
above (if the coupling is sufficiently strong), they have many channels for decay and thus 
will maintain local equilibrium through the phase transition. 

To make this argument more quantitative, consider the situation in which the final tem- 
perature is below the critical value and early times after the transition. For small amplitude 
fiuctuations of the field, long wavelength modes "see" an inverted harmonic oscillator and 
the amplitude fiuctuations begin to grow as (see below) 

{<^km-M - exp[2W{k)t] (2.3) 



W{k) = \Jfi^iT) - P (2.4) 
AT)=Am-i^?] (2.5) 

for P < fi'^{T). 

In particular this situation, modelled with the "inverted harmonic oscillators" is precisely 
the situation thoroughly and clearly studied by Guth and Pi and Weinberg and Wu pGj . 



The time scales that must be compared for the dynamics of these instabilities, are now 
the growth rate T{k) ^ \J jj? (T) — fc^ and the expansion rate H /Mpi ^ (10~^)T, if 
the expansion rate is comparable to the growth rate, then the long wavelength modes that 
are unstable and begin to correlate may be in local thermodynamic equilibrium through the 
cooling down process. Using Tc ~ /i(0)/V^, we must compare [1 — (|r)^]^ to 
Clearly for weakly coupled theories, or very near the critical temperature, the growth rate 
of the unstable modes will be much slower than the rate of cooling down, and the phase 
transition will be supercooled, similarly to a "quenching process" from a high temperature, 
disordered phase to a supercooled low temperature situation. For example, for A ~ 10^^^ the 
growth rate of long wavelength fiuctuations, is much smaller than the expansion rate even 
for a final temperature Tf ^ 0, and the long wavelength modes will be strongly supercooled. 

As mentioned previously, we do not attempt in this article to study the situation in an 
expanding gravitational background, and limit ourselves to studying the dynamics of the 
phase transition in Minkowski space by modelling the important features that are relevant 
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for the phase transition in a weakly coupled theory. Our goals here are to introduce the 
methods to study this type of phase transitions out of equilibrium, and to study the physics 
of domain formation and growth within a simplified situation. Eventually we propose to 
extend these methods to the case of an inflationary background. 

We do not envisage here to account for the cooling down process (which requires a clear 
understanding of gravitational effects, and time scales) and restrict ourselves to assuming a 
supercooled phase transition in a weakly coupled theory and propose a particular model to 
understand this situation in Minkowski space. 

III. STATISTICAL MECHANICS OUT OF EQUILIBRIUM: 

A period of rapid temperature change, may be modeled by considering a time dependent 
Hamiltonian with a time dependent mass term of the form 



H(t) 


= / d^x 
Jn 




(3.1) 


m^t) 




-t) - m}e{t) 


(3.2) 


m\ 


rp2 

c 


-1] 


(3.3) 




= - 


^] 
7-'2J 

c 


(3.4) 



with jj,^ > ; Ti > Tc ; Tf <^ T^. The introduction of the Tj ; Tf in the above mass term, is 
just a parametrization of the model, which incorporates the ingredients of a high temperature 
state at t < and a low temperature situation for t > 0. Again, the mechanism that drives 
the phase transition may either be a period of rapid inflation or a sudden coupling to a heat 
bath at a much lower temperature. The above parametrization incorporates by hand this 
"rapid supercooling" situation. 

Clearly this is a simplification, but in view of the above comments, we believe that this 
approximation is justified insofar as we are trying to understand the dynamics of the in- 
stabilities of the long-wavelength modes and the growth of correlations in weakly coupled 
theories. This assumption of a rapid "quench" may be relaxed at the expense of compli- 
cations. As it will become clear below, this approximation will allow us to obtain analytic 
results and to perform explicit calculations. Furthermore, we assume that for all times t < 
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the system is in thermal equihbrium at the initial temperature T,, thus described by the 
density matrix 

A = e-^'^' (3.5) 

Hi = H{t < 0) (3.6) 

In the Schroedinger picture, the density matrix evolves in time as 

p(t) = UmU-\t) (3.7) 

with U{t) the time evolution operator. 

An alternative, and equally valid interpretation is that we consider an initial condition in 
which the system is in thermodynamic equilibrium at temperature Tj ^ in the symmetric 
phase for t < 0, and evolved in time with a Hamiltonian with a "negative mass squared" 
that allows for broken symmetry states for t > 0. 

This interpretation in fact describes the situation studied by Guth and Pi p5[ and Wein- 
berg and Wu [Q. These authors prepare an initial gaussian state or density matrix, and 
study the time evolution of this initially prepared state with a Hamiltonian for a collection 
of "inverted harmonic oscillators" . Preparing an initial state, and evolving it with a Hamil- 
tonian of which the initial state is not an eigenstate (in the language of density matrices, the 
density matrix does not commute with the Hamiltonian) is the quantum mechanics equiva- 
lent of a "quenching process" or a "sudden approximation" . It is in this sense that we are 
thus generalizing the situation studied by the above authors. 

The expectation value of any operator is thus 

{0){t) =Tre-'^^"^U-\t)OU{t)/Tre''^^"^ (3.8) 

This expression may be written in a more illuminating form by choosing an arbitrary time 
T < for which U{T) = exp[— zTifj] then we may write exp[— = exp[— zifj(T — — 
T)] = U(T — iPi,T). Inserting in the trace U''^(T)U(T) = 1, commuting U^^(T) with pi 
and using the composition property of the evolution operator, we may write ( p.8| ) as 

{O) (t) = TrUiT - iA, t)OU{t, T)/TrU{T - if3i, T) (3.9) 
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The numerator of the above expression has a simple meaning: start at time T < 0, evolve 
to time t, insert the operator O and evolve backwards in time from t to T < 0, and along 
the negative imaginary axis from T to T — if3i. The denominator, just evolves along the 
negative imaginary axis from T to T — i(3i. The contour in the numerator may be extended 
to an arbitrary large positive time T' by inserting U{t,T')U{T' ,t) = 1 to the left of O in 
( ^I9D , thus becoming 

{0){t) = TrU{T - i(3i,T)U{T,T')U{T',t)OU{t,T)/TrU{T - i(3i,T) (3.10) 

The numerator now represents the process of evolving from T < to t, inserting the operator 
O, evolving further to T', and backwards from T' to T and down the negative imaginary 
axis to T — ij3i. This process is depicted in the contour of figure (1). Eventually we take 
T — oo ; T' oo. It is straightforward to generalize to real time correlation functions of 
Heisenberg picture operators. 

This formalism allows us also to study the general case in which both the mass and the 
coupling depend on time. The insertion of the operator O, may be achieved as usual by 
introducing currents and taking variational derivatives with respect to them. 

Because the time evolution operators have the interaction terms in them, and we would 
like to generate a perturbative expansion and Feynman diagrams, it is convenient to intro- 
duce source terms for all the time evolution operators in the above trace. Thus we are led 
to consider the following generating functional 

Z[J+, J^] = TrU{T - i(3i, T- J^)U{T, T'; J')U{T', T; J+) (3.11) 

The denominator in (|3.9| ) is simply Z[0,0,0] and may be obtained in a series expansion in 
the interaction by considering Z[0,0,J^]. By inserting a complete set of field eigenstates 
between the time evolution operators, finally, the generating functional Z[J^, , J^] may 
be written as 

Z[J+ , J- , J"] = J D^D^.D^^ J X 

^^ir'^ci^^J^] (3.12) 

with the boundary conditions <I>+(T) = $^(T-z/5i) = $ ; $+(T') = <I>-(T') = $2 ; $"(T) = 
$^(T) = $1. This may be recognized as a path integral along the contour in complex time 
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shown in Figure (1). As usual the path integrals over the quadratic forms may be evaluated 
and one obtains the final result for the partition function 

rT' 



Z[J+, J-, J^] = exp <! z dt 



exp < I 



T-il3, 



dtCintHS/dJ") I exp 1^ ^ dti ^ dt2Jc{ti)Jc{t2)Gc{ti, 
Where Jc stands for the currents on the contour as shown in figure (1), Gc are the Green's 



L7 '-2; 



(3.13) 



functions on the contour [^], and again, the spatial arguments were suppressed. 

In the limit T — oo, the contributions from the terms in which one of the currents is or 
and the other is a vanish when computing correlation functions in which the external 
legs are at finite real time, as a consequence of the Riemann-Lebesgue lemma. For this real 
time correlation functions, there is no contribution from the terms that cancel between 
numerator and denominator. Finite temperature enters through the boundary conditions on 
the Green's functions (see below). For the calculation of finite real time correlation functions, 
the generating functional simplifies to p8|,p9| 



Z[J+, J-] =expli dt 



Cint{-iS/6J^) - Cint{iS/6J~ 



X 



rT' 

'^^1 dt2Jaiti)Jti(t2)Gab(ti,t2 



(3.14) 



with a,b = +, —. 

The Green's functions that enter in the integrals along the contours in equations (|3.13 
3.14| ) are given by (see above references) 



^++(^1,^1; f2, ta) = G>{n, ti- r2, t2)Q{ti - ts) + G'<(fi, ti; ra, t2)Q{t2 - ti) 
G'~{riM] r2, t2) = G>(fi, ti, fs, ^2)6(^2 - ti) + G<(fi, ti; r2, ^2)6(^1 - ^2) 
G'+"(fi,ti;f2,t2) = -G'<(fi,ti;f2,t2) 

G-+(fi,ti;f2,t2) = -G>{fiM\r2M) = -G<{n,tr,f2,t2) 

G>(fi,ti;f2,t2) = ($(Fi,ti)$(f2,t2)) 

G< (fi, T; r2, t2) = G>(fi, T - zf3f, r2, ^2) 



(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 



The condition( |3.20| ) is recognized as the periodicity condition in imaginary time and 
is a result of considering an equilibrium situation for t < 0. The functions G^, G'^ are 
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homogeneous solutions of the quadratic form, with appropriate boundary conditions and 
will be constructed explicitly below. 

This formulation in terms of time evolution along a contour in complex time has been 
used many times in non-equilibrium statistical mechanics. To our knowledge the first to 
introduce this formulation were Schwinger and Keldysh ||^ (for an early account see 
Mills |^2|]). There are many clear articles in the literature using this techniques to study 
real time correlation functions [p^ , p9| ,p^-|37 . 



Our goal is to study the formation and growth of domains and the time evolution of 
the correlation functions. In particular, the relevant quantity of interest is the equal time 
correlation function 



s{f-t) = {<i>{f,tM6,t)) 



Sif;t) 



ik-f 



S{k;t) 



(27r)3 

S{k;t) = = HGt^t;t)) 



(3.21) 
(3.22) 
(3.23) 



where we have performed the Fourier transform in the spatial coordinates (there still is 
spatial translational and rotational invariance). Notice that at equal times, all the Green's 
functions are equal, and we may compute any of them. 

Clearly in an equilibrium situation this equal time correlation function will be time 
independent, and will only measure the static correlations. In the present case, however, 
there is a non trivial time evolution arising from the departure from equilibrium of the 
initial state. This correlation function will measure the correlations in space, and their time 
dependence. 

The function G'^{t,t') is constructed from the homogeneous solutions to the operator of 
quadratic fluctuations 



dt' 



+ k"^ + m^{t) 



ut 







(3.24) 



with IV? {t) given by (3.2). 

The boundary conditions on the homogeneous solutions are 



U^it < 0) = e^*^<('^)* 



k^ + m? 



(3.25) 
(3.26) 
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corresponding to positive frequency (particles) and negative frequency (antiparticles) 
{U^{t)]U^{t), respectively). 

The solutions are as follows: i): stable modes (/c^ > m^) 



1 / uj^jky 

2 I ^ u^^k)^ 



1 



h, = -:\l 



ii): unstable modes (/c^ < mj) 

u-k{t) = {u+k{t)Y 



W{k) = \ m) - P 



Ak = 



W{k) J 



Bk = {AkY 



With these mode functions, and the periodicity condition ( p. 20] ) we find 

i 1 



^^^^'^'^ 2u;<(A;)l-e-ft-<W 
G<{t,t')=G>{t\t) 



Utmu{t')+e-^^-<^'^U^{t)Ut{t') 



The zeroth order equal time Green's function becomes 



2iu<{k) 



coth[Au;<(A;)/2] 



for t < 0, and 



-{[1 + 2AkBk [cosh{2W{k)t) - 1]] Q{ml - P 



2uj<{k) 

[1 + 2akhk [cos(2cu>(A;)t) - 1]] 0(P - m))} coth[Ac<;<(A;)/2] 



(3.27) 
(3.28) 
(3.29) 

(3.30) 

(3.31) 
(3.32) 



(3.33) 
(3.34) 

(3.35) 

(3.36) 
(3.37) 



(3.38) 
(3.39) 



(3.40) 



(3.41) 
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for t > 0. 

The first term, the contribution of the unstable modes, reflects the growth of correlations 
because of the instabilities and will be the dominant term at long times. 



IV. ZEROTH ORDER CORRELATIONS: 

Before proceeding to study the correlations in higher orders in the coupling constant, 
it will prove to be very illuminating to understand the behavior of the equal time non- 
equilibrium correlation functions at tree-level. Because we are interested in the growth of 
correlations, we will study only the contributions of the unstable modes. 

The integral of the equal time correlation function over all wave vectors shows the familiar 
short distance divergences. From the above expression; however, it is clear that these may be 
removed by subtracting (and also multiplicatively renormalizing) this correlation function 
at t = 0. The contribution of the stable modes to the subtracted and multiplicatively 
renormalized correlation function is always bounded in time and thus uninteresting for the 
purpose of understanding the growth of the fluctuations. 

We are thus led to study only the contributions of the unstable modes to the subtracted 
and renormalized correlation function, this contribution is finite and unambiguous. 

For this purpose it is convenient to introduce the following dimensionless quantities 



iTLf rrij: 



T 2 - Tf 



rrift ; x = mf-r (4.1) 



Furthermore for the unstable modes fc^ < mj, and for initial temperatures larger than the 
critical temperature = 24/2^/ A, we can approximate coth[/3ja;<(A;)/2] ^ 2Tj/u;<(fc). Then, 



at tree-level, the contribution of the unstable modes to the subtracted structure factor (|3.23|) 
S^^\k,t) - S^^\k,0) = (l/m/)5(°)(/t,r) becomes 

S''' r) = I I (^) ^Jl + §,] [cosH2W.r) - 1] (4.2) 




- c 



ujI = k' + (4.3) 
W^ = l-K^ (4.4) 
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To obtain a better idea of the growth of correlations, it is convenient to introduce the 
scaled correlation function 

The reason for this is that the minimum of the tree level potential occurs at A$^/6mj = 1, 
and the inflexion (spinodal) point, at A$^/2mj = 1, so that T'(0,r) measures the excursion 
of the fluctuations to the spinodal point and beyond as the correlations grow in time. 

At large r (large times), the product k,'^S{k,,t) in ( ^.5| ) has a very sharp peak at Ks = 
In the region x < ^fr the integral may be done by the saddle point approximation 
and we obtain for Tj/Tc ^ the large time behavior 

P(a;,r) ^D(0,r)exp[--]-^^^^ (4.6) 



r 2 



Restoring dimensions, and recalling that the zero temperature correlation length is ^ (0) = 
l/-\/2/i, we find that for 7/ the amplitude of the fluctuation inside a "domain" ($^(t)), 
and the "size" of a domain ^d(^) grow as 

(V2«/«0))? 



2 



Ut) ^ (8^)^^(0)\/^ (4-9) 

An important time scale corresponds to the time at which the fluctuations of the 
field sample beyond the spinodal point. Roughly speaking when this happens, the insta- 



bilities should shut-off as the mean square root fluctuation of the field y($^(t)) is now 
probing the stable region. This will be seen explicitly below when we study the evolution 
non-perturbatively in the Hartree approximation and the fluctuations are incorporated self- 
consistently in the evolution equations. In zeroth order we estimate this time from the 
condition 3P(0,t) = 1, we use A = 10~^^ ; Tj/T^ = 2, as representative parameters (this 
value of the initial temperature does not have any particular physical meaning and was 
chosen only as representative). We find 

Ts ^ 10.15 (4.10) 
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or in units of the zero temperature correlation length 



t ^ 14.2^(0) (4.11) 

for other values of the parameter Tg is found from the above condition on ( [4.7|) . 
These are some of the main results of this work. 



V. PERTURBATION THEORY AND ITS DEMISE: 

The results presented in the previous section, rely on a zero-order (tree level) analysis 
of the non-equilibrium correlation function. Clearly one needs to incorporate the effects of 
the interaction. The non-equilibrium formalism introduced above lends itself to a diagram- 
matic expansion of the non-equilibrium correlation functions. We now present a one-loop 
calculation of the equal time correlation function ($^(t)$_g(i(:)). 

There are two vertices, corresponding to forward (+) and backward (— ) time propagation, 
with opposite couplings and four different propagators as given in equations ( |3.15[ - prT8[) (see 



figure (2. a)). The Feynman rules are the standard ones. The Feynman diagrams that 
contribute up to one loop to the structure factor (|3.23| ) 



are depicted in figure (2.b). Then up to one loop, and in terms of the zero order Green's 
functions, we find 

(<|.g(t)$_,Kt)) = (-zG<(t,t)) 

X ft f (Pq 



+ 2 j_Ji^ J 7^{^o>.(ii,^i) [iGldt,h)r - (G<,(t,ti)) 



X 



(2vr) 

coth[Acu<(g)/2] coth[Ac^<(A;)/2]} (5.1) 

where in the one loop integral we wrote the finite temperature Green's functions in terms 
of the zero temperature ones Gq ^. Clearly because of the complicated time dependence, 
the Dyson's series for the propagator may not be summed exactly and we only analyze here 
the one loop contribution given above. Before proceeding further with the analysis, let us 
understand the renormalizations that are necessary. 



17 



It becomes more illuminating to write the one-loop contribution explicitly in term of the 
mode functions, 

m)K^,(')) - ^"2]%)^ coth[Ac.<(fc)/2] + 



2 



X 



{u^{t)U^{t,)f -{u^{t,)U^{t))' 
colh[(5^uJ<{q)/2] coth[/3iU;<(A;)/2] (5.2) 

Clearly, in the one-loop contribution, the terms with wavevectors k correspond to the "ex- 
ternal legs" , whereas the terms with q which are integrated over, correspond to the loop line 

in figure (2.b) (because at equal time the -|— I- and terms are equal). 

First let us study the above contribution for t < 0, in which case we should recover the 
usual result. In this case both t] ti < 0, and performing the time integral with an adiabatic 
cutoff we find 

\ f d^q ( 1 \ ( 1 \^ 
' ' ' ' ' X 



'<{k) ' ^/ J 2 J (27r)3 \2u;<{q) J \2u;4k)y 

coth[PiU<{q) /2] coth[PiUJ<{k) /2] (5.3) 

In fact this is the familiar equal time Green's function up to one loop of the time in- 
dependent theory. The one loop term has temperature independent ultraviolet divergent 
contribution arising from the q-integral. Introducing an upper momentum cutoff A and an 
arbitrary renormalization scale /C we obtain 



W«.<0)=/|^(^)=^[A'-m?lu( 



•A 
X 



(5.4) 



In any case, it becomes clear that the potential divergences of the one-loop contribution, 
can be traced to the integral (again only the zero temperature contribution is divergent) 

d^q f 1 
{2nf\2^{q)^ 

For t > 0, the time integral in the one-loop correction can be split into the integral 

/°oo '^^1 ^"^^ /o "^^i- ■'-^ the first integral (from —00 to 0), {ti)U~ {ti) = 1 and the divergence 
structure is the same as that analyzed for t < In the second integral (from to i)) the 



= /7|f3(^)<(*.H-(«0 (5.5) 



divergent contribution arises solely from the stable modes as the unstable modes are cutoff 
a.t q = TUf. By analyzing the product of the mode functions, it becomes clear that the time 
dependent part {2aqbgCos[2u^{q)ti]) will yield to a finite contribution because the strong 
oscillations (for ti > 0) ensure convergence at large momenta. Thus the divergent term 
arises only from the time independent term, (a^ + 6^), in the product of mode functions. 
Finally we find the divergent term to be temperature independent and for ti > given by 



2 




ili^^W^^-A-^-JrW (5-6) 



Again in terms of an ultraviolet cutoff (A) and renormalization scale (/C), we find 

1 



hivih > 0) 



A^-(-mJ)ln(A 



(5.7) 



These divergences may be cancelled by introducing a local but time dependent countert- 
erm in the original Lagrangian density 

= 6m\t)<l>\r,t) (5.8) 
6m'{t) = / dq f {e{-t) + emiq' - m)){al + fej)} (5.9) 



On the forward and backward time contour for the non- equilibrium theory, this counterterm 
translates in the two counterterm insertions shown in figure (2.c). The introduction of these 
counterterms renders finite the one-loop contribution to all the non-equilibrium Green's 
functions as may now be easily checked. 

Having disposed of the renormalization problem, we must however address the issue 
of the instabilities. The instabilities and growth of correlations at zeroth-order had been 
analyzed before. We now realize that in the loop integral there is a contribution to the loop 
from the integration over the unstable modes which will enhance the exponential growth in 
the correlation functions. 

The maximum instability in the one-loop term is when the mode functions for both 
momenta q\k are unstable, {q^^k"^ < rrij). For the initial temperature T, > Tc, for these 
values of the momenta we use the high temperature approximation coth[/3jti;</2] ^ 2Tj/ci;<. 
It is convenient to introduce the dimensionless quantities defined in equation ( |4.1| ) and 
Q = q/rrif. Using equations ( p.3| , |3.4|) and = 24/i^/A, and the same conventions as for the 
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zeroth order structure factor given by equations ( |4.3|J4.4| , ^?TD , we obtain for the most unstable 
contribution to the one-loop correction to the structure factor S^^\k,t) = (l/mj)iS*^^^(/t, r) 



T 



n^J T2[l - (T|/T2)] Jo 
1 



dri 







X 



sinh[W,(r-ri)]{l 



1 + 



[cosh(2VrQri) - 1]} X 



1 + 



cosh[W^,(r + n)] + ( 1 - 1 cosh[W^,(r - n)] 



(5.10) 



The integral over ti may be carried out yielding a rather cumbersome result, but it becomes 
clear that this result will grow roughly as the square of the zeroth order result at large r. 
Introducing the scaled correlation function as in equation ( [4 .51 ) both for the zeroth order and 
the one-loop order 1)^^\x,t) ; V^^\x,t) respectively, in figure (3. a) we show the behavior 
for 3P^°''(0,r) (solid lines) and 3'D*-^''(0, r) (dashed lines) for the values of the parameters 
A = 10~^^ , (Tj/Tc) = 2. It is clear that eventually the one loop term becomes much larger 
than the tree level term even in the case of very weak coupling. This is a consequence of 
the instabilities and the growth of correlations that are a hallmark of the phase transition. 
Clearly the dynamics of the phase transition cannot be studied in perturbation theory. In 
fact this result in a very quantitative manner confirms the ideas that the onset of the 
transition and the time evolution of the system after the phase transition cannot be studied 
perturbatively. 



VI. BEYOND PERTURBATION THEORY: HARTREE APPROXIMATION 

It became clear from the analysis of the previous section that perturbation theory is in- 
adequate to describe the non-equilibrium dynamics of the phase transition, precisely because 
of the instabilities and the growth of correlations. This growth is manifest in the Green's 
functions that enter in any perturbative expansion thus invalidating any perturbative ap- 
proach. Higher order corrections will have terms that grow exponentially and faster than 
the previous term in the expansion. And even for very weakly coupled theories, the higher 
order corrections eventually become of the same order as the lower order terms. 

As the correlations and fluctuations grow, field configurations start sampling the stable 
region beyond the spinodal point. This will result in a slow down in the growth of cor- 
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relations, and eventually the unstable growth will shut-off. When this happens, the state 
may be described by correlated domains with equal probabibility for both phases inside the 
domains. The expectation value of the field in this configuration will be zero, but inside each 
domain, the field will acquire a value very close to the value in equilibrium at the minimum 
of the effective potential. The size of the domain in this picture will depend on the time 
during which correlations had grown enough so that fluctuations start sampling beyond the 
spinodal point. 

Since this physical picture may not be studied within perturbation theory, we now intro- 
duce a non-perturbative method based on a self-consistent Hartree approximation []38|-^ . 



The self-consistent Hartree approximation is implemented as follows: in the initial La- 
grangian write 

A$4(f,t) = ^($^(f,t))<|.2(f,t) + (^$^(r,t) - ^($2(f,t))<|.2(f,t)^ (6.1) 

the first term is absorbed in a shift of the mass term 

m'{t)^m'{t) + ^{^\t)) 

(where we used spatial translational invariance). The second term in (|6.1|) is taken as an in- 
teraction with the term ($^(t))$^(r, t) as a "mass counterterm" . The Hartree approximation 
consists of requiring that the one loop correction to the two point Green's functions must 
be cancelled by the "mass counterterm" . This leads to the self consistent set of equations 

{^\t)) = / ^ {-^GUt,t)) = J ^.i^^^tmkit) coth[A^<(A:)/2] (6.2) 



W± = (6.3) 



Before proceeding any further, we must address the fact that the composite operator 
($^(7^, t)) needs one subtraction and multiplicative renormalization. As usual the subtraction 
is absorbed in a renormalization of the bare mass, and the multiplicative renormalization 
into a renormalization of the coupling constant. We must also point out that the Hartree 
approximation is uncontrolled in this scalar theory; it becomes equivalent to the large-N 
limit in theories in which the field is in the vector representation of 0(N) (see for example 

i)- 
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At this stage our justification for using this approximation is based on the fact that it 
provides a non-perturbative framework to sum an infinite series of Feynman diagrams of the 
cactus type 

In principle one may improve on this approximation by using the Hartree propagators 
in a loop expansion. The cactus-type diagrams will still be cancelled by the counterterms 
(Hartree condition), but other diagrams with loops (for example diagrams with multiparticle 
thresholds) may be computed by using the Hartree propagators on the lines. This approach 
will have the advantage that the Hartree propagators will only be unstable for a finite time 
t < ts- It is not presently clear to these authors, however, what if any, would be the 
expansion parameter in this case. 

It is clear that for t < there is a self-consistent solution to the Hartree equations with 
equation (|6.2| ) and 

($^(t)) = ($2(0-)) 

= explTtoo^ik)] (6.4) 



u;i{k) = k' + mt + - + ($2(0-)) = k' + ml^ 



(6.5) 



where the composite operator has been absorbed in a renormalization of the initial mass, 
which is now parametrized as m^^ = fij^liT^ /T^) — 1]. For t > we subtract the compos- 
ite operator at t = absorbing the subtraction into a renormalization of m'j which we now 
parametrize as m'j = — (Tj/T^)]. We should point out that this choice of parametriza- 
tion only represents a choice of the bare parameters, which can always be chosen to satisfy 
this condition. The logarithmic multiplicative divergence of the composite operator will be 
absorbed in a coupling constant renormalization consistent with the Hartree approximation 



41| , however, for the purpose of understanding the dynamics of growth of instabilities as- 



sociated with the long-wavelength fluctuations, we will not need to specify this procedure. 
After this subtraction procedure, the Hartree equations read 

[($2(t)) - ($2(0))] = J ^.l^^.mmuit) - 1] coth[Au;<(fc)/2] (6.6) 



^ + P + m^(t) + :^(($2(t))-($2(0)) 



Ut{t) = (6.7) 
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2 



T- 

^ - 1 

2^2 



1 — ^- 

2^2 



e(t) (6.8) 



with Tj > Tc and Ty -C Tc. With the self-consistent solution and boundary condition for 
t < 

< 0)) - ($'(0))] = (6.9) 
U^{t<0) = exp [TiuJ< {k)t] (6.10) 



u<{k) = ^ + mj^ (6.11) 

This set of Hartree equations is extremely complicated to be solved exactly. However 
it has the correct physics in it. Consider the equations for t > 0, at very early times, 
when (the renormalized) (^^(t)) — ($^(0)) ~ the mode functions are the same as in the 
zeroth order approximation, and the unstable modes grow exponentially. By computing 
the expression ( |6.6| ) self-consistently with these zero-order unstable modes, we see that the 
fluctuation operator begins to grow exponentially. 

As ((^^(t)) — ($2(0))) grows larger, its contribution to the Hartree equation tends to 
balance the negative mass term, thus weakening the unstabilities, so that only longer wave- 
lengths can become unstable. Even for very weak coupling constants, the exponentially 
growing modes make the Hartree term in the equation of motion for the mode functions 
become large and compensate for the negative mass term. Thus when 



(($^(t)) - ($^(0)) 



the instabilities shut-off, this equality determines the "spinodal time" tg- The modes will 
still continue to grow further after this point because the time derivatives are fairly (expo- 
nentially) large, but eventually the growth will slow-down when fluctuations sample deep 
inside the stable region. 

After the subtraction, and multiplicative renormalization (absorbed in a coupling con- 
stant renormalization) , the composite operator is finite. The stable mode functions will 
make a perturbative contribution to the fluctuation which will be always bounded in time. 
The most important contribution will be that of the unstable modes. These will grow expo- 
nentially at early times and their effect will dominate the dynamics of growth and formation 
of correlated domains. The full set of Hartree equations is extremely difficult to solve, even 
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numerically, so we will restrict ourselves to account only for the unstable modes. From the 
above discussion it should be clear that these are the only relevant modes for the dynam- 
ics of formation and growth of domains, whereas the stable modes, will always contribute 
perturbatively for weak coupling after renormalization. 

Introducing the dimensionless ratios ( |4.1j ) in terms of rrif^R ; mj^/j, (all momenta are 
now expressed in units of Jtj-/,/?), dividing ( |6.7|) by rrij j^, using the high temperature approx- 
imation coth[j3iUJ^{k) /2] ^ 2Ti/uj^{k) for the unstable modes, and expressing the critical 
temperature as T"^ = 24/iij/Ai?, the set of Hartree equations ( |6.6| , [HT^ ) become the following 
integro-differential equation for the mode functions for t > 



^ + 9 dp- 



LI 



[u;{t)u^{t)-i\ 



U^{t) = 



with 



U^{t<0) = exp[T^u;<{q)t] 
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Tf] 



R 



T, 



(6.12) 

(6.13) 
(6.14) 

(6.15) 
(6.16) 



47r2 [7^2 

The effective coupling (|6.16|) reflects the enhancement of quantum fluctuations by high tem- 
perature effects; for Tf/Tc ~ 0, and for couplings as weak as Xr ~ 10^^^, g ~ 10^'^(Tj/Tc). 

The equations ( |6.12| ) may now be integrated numerically for the mode functions; once we 
find these, we can then compute the contribution of the unstable modes to the subtracted 
correlation function equivalent to ( [4.5|) 

Xr 



3I)(^^)(x,r) 

In figure (4) we show 
Xr 



(<l>(f,t)<l>(0,t))-(<l>(f,0)<l>(0,0)) 



9 I dp 
lo 



p^ 



sm{px) -. 
(px) 



(6.17) 
(6.18) 



((<|.^(r)) - ($^(0))) = 3(P^^(0, r) - (0, 0)) 



(solid line) and also for comparison, its zeroth-order counterpart 3(P*-°''(0, r) — P*^°)(0,0)) 
(dashed line) for \r = 10"^^ , Ti/Tc = 2. (Again, this value of the initial temperature does 
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not have any particular physical significance and was chosen as a representative) . We clearly 
see what we expected; whereas the zeroth order correlation grows indefinitely, the Hartree 
correlation function is bounded in time and oscillatory. At r 10.52 , 3(P*^^^)(0, r) — 
X)(^^)(0,r)) = 1, fluctuations are sampling field configurations near the classical spinodal, 
fiuctuations continue to grow, however, because the derivatives are still fairly large. However, 
after this time, the modes begin to probe the stable region in which there is no exponential 
growth. At this point (($^(t)) — $^(0))), becomes small again because of the small 
coupling g ^ 10~^, and the correction term becomes small. When it becomes smaller than 
one, the instabilities set in again, modes begin to grow and the process repeats. This gives 
rise to an oscillatory behavior around (($^(t)) — $^(0)) = 1 as shown in figure (4). 

f,R 

In figures (5.a-d), we show the structure factors as a function of x ior t — 6, 8, 10, 12, 
both for zero-order (tree level) X>(°^ (dashed fines) and Hartree (sofid fines). These 

correlation functions clearly show that correlations grow in amplitude and that the size of 
the region in which the fields are correlated increases with time. Clearly this region may 
be interpreted as a "domain", inside which the fields have strong correlations, and outside 
which the fields are uncorrelated. 

We see that up to the spinodal time ^ 10.52 at which 2^^(($^(rs)) - $^(0))) = 1, 
the zeroth order correlation 3P*^°^(0, r) is very close to the Hartree result. In fact at r^, 
the difference is less than 15% In particular for these values of the coupling and initial 
temperature, the zeroth order correlation function leads to ts ~ 10.15, and we may use the 
zeroth order correlations to provide an analytic estimate for tg, as well as the form of the 
correlation functions and the size of the domains. 

The fact that the zeroth-ordcr correlation remains very close to the Hartrcc-corrected 
correlations up to times comparable to the spinodal is clearly a consequence of the very 
small coupling. 

To illustrate this fact, we show in figures (6,7) the same correlation functions but for 
A — 0.01, Tj/Tc — 2. Clearly the stronger coupling makes the growth of domains much faster 
and the departure from tree-level correlations more dramatic. Thus, it becomes clear that 
for strong couplings, domains will form very rapidly and only grow to sizes of the order of the 
zero temperature correlation length. The phase transition will occur very rapidly, and clearly 
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our initial assumption of a rapid supercooling will be unjustified. This situation for strong 
couplings, of domains forming very rapidly to sizes of the order of the zero temperature 
correlation length, is the picture presented by Mazenko and collaborators However, 
for very weak couplings (consistent with the bounds from density fluctuations), our results 
indicate that the phase transition will proceed very slowly, domains will grow for a long time 
and become fairly large, with a typical size several times the zero temperature correlation 
length. In a sense, this is a self consistent check of our initial assumptions on a rapid 
supercooling in the case of weak couplings. 

Thus, as we argued above, for very weak coupling, we may use the tree level result to 
give an approximate bound to the correlation functions up to times close to the spinodal 
time using the result given by equation (|4.7|), for Tf ^ 0. Thus, we conclude that for large 
times, and very weakly coupled theories {Xr < 10^^^) and for initial temperatures of the 
order of the critical temperature, the size of the domains ^oit) will grow typically in time 
as 



with ^(0) the zero temperature correlation length. The maximum size of a domain is 
approximately determined by the time at which fluctuations begin probing the stable region, 
this is the spinodal time tg and the maximum size of the domains is approximately ^nits)- 

An estimate for the spinodal time, is obtained from equation ( |4.7] ) by the condition 
3V{ts) = 1, then for weakly coupled theories and Tj ^ 0, we obtain 



To = —;=^ — In 

v^^(O) 



47r3 



/ ( Tj 



(6.20) 



It is remarkable that the domain size scales as ^d(^) ~ just like in classical theories 
of spinodal decomposition, when the order parameter is not conserved, as is the case in 
the scalar relativistic field theory under consideration, but certainly for completely different 
reasons. At the tree level, we can identify this scaling behavior as arising from the relativistic 
dispersion relation, and second order time derivatives in the equations of motion, a situation 
very different from the classical description of the Allen- Cahn- Lifshitz theory of 



spinodal decomposition based on a time-dependent Landau Ginzburg model. 
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Beyond Hartree: 

The Hartree approximation, keeping only the unstable modes in the self- consistent equa- 
tion, clearly cannot be accurate for times beyond the spinodal time. When the oscillations in 
the Hartree solution begin, the field fluctuations are probing the stable region. This should 
correspond to the onset of the "reheating" epoch, in which dissipative effects become impor- 
tant for processes of particle and entropy production. Clearly the Hartree approximation 
ignores all dissipative processes, as may be understood from the fact that this approximation 
sums the cactus type diagrams for which there are no multiparticle thresholds. Furthermore, 
in this region, the contribution of the stable modes to the Hartree equation becomes impor- 
tant for the subsequent evolution beyond the spinodal point and clearly will contribute to 
the "reheating" process. A possible approach to incorporate the contribution of the stable 



modes may be that explored by Avan and de Vega [^] in terms of the effective action for 
the composite operator. 

Thus, although the Hartree approximation may give a fairly accurate picture of the 
process of domain formation and growth, one must go beyond this approximation at times 
later than the spinodal time, to incorporate dissipative effects and to study the "reheating" 
period. Clearly, one must also attempt to study the possibility of "percolation of domains" . 
Furthermore, the Hartree approximation is essentially a Gaussian approximation, as the 
wave-functional (or in this case the functional density matrix) is Gaussian with kernels that 
are obtained self- consistently. The wave-functional must include non-gaussian correlations 
that will account for the corrections to the Hartree approximation and will be important to 
obtain the long time behavior for t > tg- 



VII. CONCLUSIONS AND LOOKING AHEAD: 



The motivations of this work were twofold. First we pointed out that the dynamics of 
typical phase transitions in weakly coupled theories must be studied away from thermody- 
namic equilibrium, and introduced the methods and techniques of non-equilibrium quantum 
statistical mechanics to study this situation. 

Second we studied both analytically and numerically the case of a strongly supercooled 
phase transition in which the system initially in thermal equilibrium at an initial temperature 
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larger than the critical is cooled down to temperatures well below the transition temperature. 
The motivation here was to model a period of rapid inflation in a weakly coupled theory and 
to study the formation and growth of correlated domains. We indicated that the dynamics of 
the phase transition cannot be studied within perturbation theory because of the instabilities 
that drive the process of domain formation and growth, that is spinodal decomposition. 

We used a non-perturbative self-consistent Hartree approximation to study the time 
evolution of domain growth, reflected in the equal time two point correlation function 



with ^(0) the zero temperature correlation length. The domains, however, will grow up to a 
maximum time at which the fluctuations begin sampling the stable region. This maximum 
"spinodal time" is approximately given for weakly coupled theories by 



When the self-couplings are strong, the phase transition proceeds rapidly, and domains 
will not have time to grow substantially, and their sizes will be of the order of the zero 
temperature correlation length. 

In principle the "sudden approximation" (quenching) may be relaxed at the expense of 
complications, however the formalism presented in this paper is completely general, once the 
initial state is specifled and the boundary conditions for the mode functions are understood, 
the time evolution of correlation functions may be studied numerically. 

Clearly, the next step is to study the dynamics of the phase transition in FRW cosmolo- 
gies. In this case, there are several physical effects that will play a very important role in 
the dynamics. In particular, the red-shift of physical wave vectors, will tend to enhance the 
instabilities, because more wave vectors arc entering the unstable region as time evolves. On 
the other hand, the presence of a horizon, and the "friction" term in the Heisenberg equa- 
tions of motion, will prevent domains from growing bigger than the horizon size. Thus there 




immm . exp[y2t/e(o)] , , sm(v^./e.(t)) 




(7.1) 



(7.2) 




(7.3) 
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seems to be a competition between the different time scales tliat must be studied carefully 
to obtain any meaningful conclusion about formation and growth of correlated domains in 
FRW cosmologies. 
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Figure Captions: 

Figure 1: Contour in complex time plane to evaluate the generating functional for non- 
equilibrium Green's functions. 

Figure 2. a: Two vertices and four propagators generate the Feynman diagramatic ex- 
pansion for non-equilibrium Green's functions. 

Figure 2.b: Diagrams that contribute up to one loop to {^^(t)^_^(t')). 

Figure 2.c: Two mass counterterms. 

Figure 3: Zero and one loop contribution to the structure factor. The sohd hne repre- 
sents 3V^'^\0,t), the dashed hne represents 3V^^\0,t). 

Figure 4: Hartree (sohd line) and zero order (dashed hne) results for ^^((^^ ('?')) — 
($2(0))) = V{0, t), for A = 10-'^, 1 = 2. 

Figure 5. a: Scaled correlation functions for r = 6, as function of x, 'D^^^\x,t) (solid 
line), and V^^\x,t) (dashed line). A = 10-l^ § = 2- 

Figure 5.b: Scaled correlation functions for r = 8, as function of a;, T)^^^\x,t) (solid 
line), and V^°\x,t) (dashed line). A = 10"^^ § = 2. 

Figure 5.c: Scaled correlation functions for r = 10, as function of x, V^^^\x,t) (solid 
line), and V^^\x,t) (dashed line). A = 10-^^^ § = 2. 

Figure 5.d: Scaled correlation functions for r = 12, as function of x, V^^^\x^t) (solid 
line), and V^^\x,t) (dashed line). A = 10"^^ ^ = 2- 

Figure 6: Hartree (sohd line) and zero order (dashed hne) results for ^^((^^(t)) — 
{^'^{t))) = 3P(0,t), for A = 0.01, g = 2. 

Figure 7: Scaled correlation functions for r — 4.15, as function of x, V'^^^\x,t) (solid 
line), and V^^\x,t) (dashed line). A = 0.01, J = 2. 
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